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Using QCD sum rules with nonlocal condensates, we show that the distribution amplitude of 
the longitudinally polarized p-meson may have a shorttailed platykurtic profile in close analogy to 
our recently proposed platykurtic distribution amplitude for the pion. Such a chimera distribution 
de facto amalgamates the broad unimodal profile of the distribution amplitude, obtained with a 
Dyson-Schwinger equations-based computational scheme, with the suppressed tails characterizing 
the bimodal distribution amplitudes derived from QCD sum rules with nonlocal condensates. We 
argue that pattern formation, emerging from the collective synchronization of coupled oscillators, 
can provide a single theoretical scaffolding to study unimodal and bimodal distribution amplitudes 
of light mesons without recourse to particular computational schemes and the reasons for them. 
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I. INTRODUCTION 

Many theoretical models exist to describe the valence 
parton distribution amplitude (DA) of the 7r, the p (both 
longitudinally and transversally polarized), and other 
mesons. In particular the pion, the lightest bound state 
within Quantum Chromodynamics (QCD), provides a 
suitable “laboratory” for testing new ideas and tech¬ 
niques to catch the main ingredients of the underlying 
quark-gluon dynamics entering exclusive QCD processes 
(see pj-lg for reviews). In this paper we will present an 
amplification of the recent analysis by one of us in |J], 
continued in Q, for the pion and consider its extension 
to the longitudinally polarized p meson (p||). The new 
mode of thought in [J] is based on the ubiquitous phe¬ 
nomenon of synchronization (Sync for short) in complex 
systems and we will expand the status of this subject 
towards a deeper understanding of the meson DAs. In 
this way, we will redetermine the p\\ DA using QCD sum 
rules with nonlocal condensates (NLC) within the ap¬ 
proach developed in . 

Our primary findings to be discussed later can be 
summarized as follows: First, we provide more details 
about the structure of the Sync-inspired shorttailed (i.e., 
endpoint-suppressed) platykurtic (pk) 7r DA, proposed 
in [4|, and quantify the uncertainties of its expansion 
coefficients. This DA is a kind of chimera 1 state be¬ 
cause it mimics within the NLC-based approach char¬ 
acteristics pertaining to the dynamical chiral symmetry 
breaking (DCSB) described in terms of Dyson-Schwinger 
equations (DSE) [iil. ITo|. In particular it conserves NLC- 
generated endpoint suppression of the tt DA combining 
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1 An imaginary creature in Greek mythology made up of different 
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it with a broad downward concave shape in the central 
region. Second, using similar mathematical techniques, 
we obtain within the reliability range of the NLC ap¬ 
proach a regime of DAs for the p\\ meson characterized 
by a shorttailed platykurtic profile. Third, we employ 
statistical measures, like the kurtosis, to classify meson 
DAs with respect to their peakedness relative to the tail 
flatness and heaviness. 

The rest of the paper is organized as follows. The next 
two sections (Sec. Inland Sec. IIIII) discuss the theoretical 
basis for the description of meson DAs within QCD. Sec¬ 
tion [IV] sketches the derivation of the 7r and p || DAs from 
QCD sum rules with NLCs. We will then proceed to in¬ 
vestigate how meson DAs can be analyzed in terms of 
synchronization concepts (Sec. |V]). Synthetic meson DAs 
will be considered in Sec. IVI1 where also the important 
chimera DAs with a shorttailed platykurtic profile will 
be presented. Finally, Sec. IVIII will be reserved for the 
summary of our main results and conclusions. 


II. MESON DAs IN QCD 

Let us consider the pion DA, starting with its defi¬ 
nition. 2 Applying collinear factorization in QCD, the 
7r DA of leading-twist two, ip^\x,p 2 ), encodes the dis¬ 
tribution of the longitudinal momentum of the pion be¬ 
tween its two valence constituents: the quark and the 
antiquark, with corresponding longitudinal-momentum 
fractions x q = x = ( k° + fc 3 )/(P° + P 3 ) = k + /P + and 
x q = 1 — x = x, respectively. Its momentum-scale de¬ 
pendence stems from the renormalization of the current 


2 The exposition to follow relies upon the review in Q to which 

we refer for details and the original references. 
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operator in 


( 0 |g(^) 7 ^ 75 [^ ; 0 ]g( 0 )| 7 r(P))| 2 2 =0 = if n P p J dxe lx ^ z ' p) 

xtpM (x,n 2 ) , (1) 

where the gauge link [ 2 , 0 ] = 

Pexp[— ig f~ dWfj,A^(w)t a \ = 1 is set equal to unity on 
account of the lightcone gauge A - n = A + = 0 (n 2 = 0). 
The pion DA is related to the Bethe-Salpeter wave 
function ip^(x, k±) by integrating over the transverse 
parton momentum k±, i.e., 



Because the dependence on the momentum scale of any 
meson DA is controlled bv the Efremov-Radyushkin 
Brodsky-Lepage (ERBL) [ll], |l2| evolution equation, 
each meson DA can be expressed in terms of the one- 
loop eigenfunctions of this equation, i/) n (x) = 62(1 — 
x)Cn ' {2x — 1), with the asymptotic (asy) DA being 

given by ip^ sy (x) = 6x(l — x) = 6xx, and Cn^ 2 \2x — 1 ) 
denoting the Gegenbauer polynomials of order 3/2 within 
the complete and orthonormal basis on 1 £ [ 0 , 1 ] with 
respect to the weight xx. Thus, one has the (scale- 
dependent) conformal expansion 

OO 

</4 2 ) (x,/r 2 )= a n (y 2 )^ n {x) (3) 

n=0,2,4,... 

in terms of the nonperturbative coefficients a n (y 2 ). By 
virtue of the normalization condition dxyffl (x, /i 2 ) = 
1 , ao = 1 at any scale y 2 . 


III. GEGENBAUER DA REPRESENTATIONS 

In our approach [9, based on QCD sum rules with 
nonlocal condensates |13l - ll8j , we calculated the moments 

(£, N )* = [ dx(2x-l) N ip {2 \x,y 2 ) (4) 

do 

up to N = 10 together with their intrinsic theoretical un¬ 
certainties at the typical hadronic scale y 2 ss 1.35 GeV 2 
0. Detailed estimates of the moment uncertainties can 
be found in OSH (see also Table!]}). This scale repre¬ 
sents the average value of the Borel parameter M 2 in the 
stability window of the sum rule, notably, M 2 ~ [m 2 , so], 
where Sq is the continuum threshold and m p is the phys¬ 
ical mass of the p meson. Note that the moments (£, n )tv 
coincide by construction with the central moments 

Vn[<Ptt] = [ dx(x - vlv*]) 1 *<Pn(x) = 2~ N (£ N ) n (5) 

J 0 
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FIG. 1: (color online). The upper panel shows some char¬ 
acteristic pion DAs. The larger shaded (green) area displays 
the region of the bimodal Bakulev-Mikhailov-Stefanis (BMS) 
family of DAs derived in @] from QCD sum rules with NLCs 
using the nonlocality parameter A 2 = 0.4 GeV 2 , defind in Sec. 
m The family of the 7r platykurtic DAs is shown as a narrow 
shaded strip in red color, obtained with A 2 = 0.45 GeV 2 at 
the edge of the NLC regime. The solid lines within both 
regions denote, respectively, the BMS model from JOj and 
the pk DA discussed in the text and in 0,0- The broken 
lines show the unimodal DSE-DB DA (dashed), the DSE-RL 
DA (dashed-dotted) (both from 0), and the asymptotic DA 
(dashed-dotted-dotted). The lower panel illustrates various 
PH DAs. The larger blue shaded area contains the family of 
DAs obtained in Q using QCD sum rules with NLCs and 
A 2 = 0.4 GeV 2 , while the narrower strip in its interior in¬ 
dicates the platykurtic regime of these DAs. The solid lines 
within each band denote, respectively, the bimodal DA from 
[§!] (lowest (blue) solid line) and the platykurtic DA (upper 
solid line) derived in this work. The lower dashed line repre¬ 
sents the DA obtained from the DSE approach 0 , whereas 
the dashed-dotted-dotted line displays again the asymptotic 
DA. All DAs in both panels refer to the scale p 2 = 4 GeV 2 af¬ 
ter two-loop ERBL evolution, provided the initial proprietary 
scale was lower than this. 


of ip„(x), where 



dxxtpir(x) 


1 

2 


( 6 ) 


is the mean of the DA. Using standard techniques (see, 
for example, [l, 0 ), we extracted from the moments 
the corresponding conformal coefficients a n entering Eq. 
©• These quantities contain nonperturbative informa¬ 
tion and implicitly depend on the finite virtuality of the 
vacuum quarks, the latter expressed by means of the 
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nonlocality parameter A 2 = (qigG^a^q)/2(q(0)q(0)} w 
[0.35 — 0.45] GeV 2 . It was found that the first two coef¬ 
ficients (see Appendix C in 0 ) 

° 2 = ^(5(e 2 )-l) , (7) 

« = ¥0« 4 >-I<g + s) < 8 > 

dominate. Their values have been calculated with con¬ 
trolled accuracy in [§j. They read 02 (g 2 ~ 1 GeV 2 ) = 
0.20 and a^g 2 ~ 1 GeV 2 ) = —0.14, whereas the next 
higher coefficients were computed in the same work 
as well and found to be much smaller, viz., a 6 « 
Q-2/3; Os ~ 02/4; a 10 « 02/5, but bearing large un¬ 
certainties. Their inclusion can add refinements to the 
method, as we have discussed in j22|. This apparent hier¬ 
archy, with each subsequent coefficient becoming smaller 
with the order of the conformal expansion, is not fol¬ 
lowing from general principles; it is an inherent ele¬ 
ment of our specific approach. Indeed, one can even 
have an inverse ordering of the coefficients a n — see 
[23! for such DAs. In the final analysis, the pion DA 
at the scale g 2 > 1 GeV 2 can be written in the form 
(£ = 2x — 1 = x — x) 


<P* MS ( X ) = 6xx [l + a 2 C2 3/2) (£) + a 4 cf /2) (£)l , (9) 


where the label means Bakulev, Mikhailov, Stefanis [fj. 
This simple model probably offers a biased picture of the 
pion structure but its chief predictions are in good agree¬ 
ment with measurements and various lattice simulations, 
as detailed in a series of papers HH3, HiiHIH Idil [2!)l | . The 
reason is that physical observables, like the pion-photon 
transition or the pion’s electromagnetic form factor are 
given in terms of integrals of the DAs with smooth coeffi¬ 
cient functions. Because the leading-order anomalous di¬ 
mensions of the involved operators in the matrix elements 
between the meson state and the vacuum are positive (ex¬ 
cept 70 = 0), the higher coefficients are logarithmically 
suppressed at large scales, so that only ipo(x) = p^ y (x) 
survives. This is particularly visible in the inverse mo¬ 
ment of the pion DA, cf. Eq. m, in which the oscil¬ 
lating terms are washed out and strongly suppressed as 
the momentum increases. The profiles of the meson DAs, 
considered in this work, are shown in the upper (71) and 
the lower (py) panel of Fig. [TJ respectively. 

To leading logarithmic accuracy, the conformal coef¬ 
ficients a n are multiplicatively renormalizable and the 
anomalous dimensions are known in closed form, see, 
for example, [3J. At the next-to-leading-order (NLO) 
level, the momentum-scale dependence of the conformal 
coefficients (or moments) is more complicated owin g to 
the mixing of the operators under renormalization [301 
[32} . The diagonal elements of the corresponding two- 
loop anomalous-dimension matrix coincide with the fla¬ 
vor nonsinglet anomalous dimensions known from deeply 
inelastic scattering. They have been computed in [331 ]. 


corrected in [34[, and verified in [35[. The NLO ERBL 
kernel was calculated in [3614381 ] . while the analytic ex¬ 
pressions for the mixing coefficients to obtain the corre¬ 
sponding eigenfunctions were given in [31, -32}. Bear in 
mind that the next-to-leading-order corrections under a 
change of scale using a running coupling appear as a two- 
loop contribution of the eigenvalues and as an a s correc¬ 
tion to the eigenfunctions. A detailed exposition of the 
ERBL evolution of the pion DA at the two-loop level, as 
used in the present work, is provided in Appendix D in 


It is convenient to employ another Gegenbauer repre¬ 
sentation (“Gegenbauer-a”) proposed in I S3, notably, 


= f{{a,a%,...,a?},x) = V’o a) ( x ) 


(<*)/ 


]T a^{g 2 )^\x), 

i—2,4,... 


( 10 ) 


where 


4 “\ x ) = N a (xxr~C^(2x - 1) (11) 

and N a = 1/B(a + 1/2, a + 1/2), cn_ = a — 1/2, where 
B(x,y) is the Euler beta function. The Gegenbauer 
polynomials Cn*\2x — 1) form an orthonormal set over 
x £ [0,1] with respect to the weight [xx) a -. The dif¬ 
ference to the conformal expansion in Eq. © is that 
the order of the Gegenbauer polynomials is not a pri¬ 
ori fixed to the value 3/2, but is allowed to vary in or¬ 
der to accelerate the reconstruction procedure of meson 
DAs on 1 £ [0,1]. However, expansion (10) is not di¬ 
rectly amenable to ERBL evolution because the functions 
i/>n(x) are not its eigenfunctions. To evolve (f^k\x, g 2 ), 
expressed via CED, to another scale Q 2 > p 2 , one has 
to project it first onto the conformal basis {ip n {x)} and 
then determine a_ and a“ at the new scale. The authors 
of the works in Refs. iE3 find that it is sufficient to 
include only one coefficient in this expansion, namely, af, 
so that Eq. m reduces to 

(x,g 2 ) = N a ( xx) a -[l + a^C^ix - x)\. (12) 

Below, we will use for our analysis both representations 
in parallel and present the results in the form (02,04) 
and (N a , a_, )■ It is worth bearing in mind that broad 

DAs of the form p v {x) ~ ( xx) a ~ are not well represented 
by Eq. © which only employs the first two conformal 
coefficients. To approximate such DAs with sufficient 
accuracy, one would have to include a large number of 
coefficients of order 50 or more, see, for example, [31 ]. 


IV. 7r AND pi, DAs FROM NLC QCD SUM 
RULES 

In this section, we discuss the derivation of the py DA 
using QCD sum rules with nonlocal condensates within 
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FIG. 2: (color online). Computed locations of various py and n DAs at p 2 = 4 GeV 2 in the plane spanned by the conformal 
coefficients a 2 and C 14 (left panel) and corresponding results in the plane (a-jdf) (right panel). The slanted (blue) rectangle 
further to the left shows the region of { 02 , 04 } values obtained with NLC-QCD sum rules and A 2 = 0.4 GeV 2 for the py DA in 
[§), with the smaller (red) rectangle in its interior denoting the platykurtic regime determined in this work. The other small 
(light-green) rectangle next to it contains the platykurtic region calculated for the pion DA in this work using A 2 = 0.45 GeV 2 
at the edge of the NLC regime. The larger slanted (green) rectangle further to the right contains the original values of { 02 , 04 } 
for the pion DA obtained with A 2 = 0.4 GeV 2 in | 6 j. The symbols ®, ©, and ® show examples of synthetic DAs according 
to Eq. (1261) for a = 0.25, a = 0.5, and a = 0.75, respectively. Corresponding results for the py DA at 4 GeV 2 in the (a_,af) 
plane are displayed in the right panel with the following designations: The stretched (blue) strip shows the region of (a_, aj) 
values calculated with NLC-QCD sum rules and A 2 = 0.4 GeV 2 in this work. The two smaller slanted rectangles illustrate the 
platykurtic regime for the py DA (upper red rectangle) and for the 7 r DA (lower light-green rectangle); they correspond to the 
analogous graphs in the left panel. The other symbols in both panels denote the following DAs: ★ bimodal py DA from 
NLC-QCD sum rules @]; ♦ bimodal py DA from NLC-QCD sum rules [3; © — platykurtic py DA (this work); T - py 
DA from lightfront model [3!| ; • - py DA from instanton model |40| ; ■ - py DA from QCD sum rules [4d| : □ - py DA 

from AdS/QCD Q A — tt AdS/QCD 0; 0 - py DSE DA 0; ♦ - asymptotic DA; X (both panels) — BMS tt DA 0; 

(both panels) — platykurtic n DA; 0]; ▲ (both panels) — DSE-DB tv DA Jt|; v (both panels) - DSE-RL 7r DA ft|. The 
slanted rectangles, bounded by dashed lines, in both panels display the regions determined in this work for the 7r DA using 
NLC QCD sum rules with A 2 = 0.45 GeV 2 . When the proprietary scale of the studied DAs was different from 2 GeV, NLO 
evolution was applied to determine their coefficients at this scale. The numerical values of all coefficients are given in Table [1] 


the scheme developed in . The main conceptual idea 
is to apply the sum-rule method in combination with 
vacuum averages of nonlocal operators fTTH f]~8l | . Following 
this rationale, one can determine the moments (£, N ) of a 
meson DA using a sum rule with nonlocal condensates. 
In addition, due to the absence of endpoint singularities 
in the NLC approach, one is able to calculate the inverse 
moment 

(a; -1 )^ = / dxx~ 1 tp n (x) (13) 

J 0 

via an independent sum rule at the same low renormal¬ 
ization scale /i > 1 GeV. This moment is an integral 
characteristic of the pion DA and encodes information on 
the maximum possible weight of the higher-order confor¬ 
mal coefficients. Moreover, it is particularly relevant for 
phenomenological applications because 

{x~X = 3(1 + 02+04 ...) = -A-Q 2 F^ 0 (Q 2 ), (14) 

where F^^ 0 (Q 2 ) is the leading-order (LO) expression of 
the pion-photon transition form factor, which has been 
measured in several experiments from a few GeV 2 up 
to 40 GeV 2 [i3 - [47j . Here the ellipsis represents correc¬ 
tions due to higher eigenfunctions and evolution of the 


Gegenbauer coefficients a n to the considered scale Q 2 is 
assumed (see [g] for more details). Thus, experimental 
evidence can be used to validate or reject particular pion 
DAs. For instance, we know (see, e.g., [0) that all exist¬ 
ing data demand (ir -1 )^ > 3, implying that the 7 r DA has 
to be broader than <p asy at accessible momentum values. 

To obtain the 7r DA, we employ the following QCD 
sum rule (m„ = 0, f w = 0.137 GeV) 


f%ip n {x)e m ” /M2 + v+bO) e m ^i /M ' 


s o 

J Ppert (s, x) e~ s/M2 ds + $ % (x, M 2 ), (15) 


whereas for the p\\ DA we have 



So 

= JPpert ( 5 , x) e~ s/M2 ds + $p(z, M 2 ) (16) 

0 


with m p = 0.775 GeV and m p > = 1.496 GeV. The ef¬ 
fective Ai-meson state with the decay constant f= 
0.227 GeV and the mass toa x = 1.2712 GeV comprises 
the n' and the ai mesons and is described by the DA 
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ipA 1 ( x ). The perturbative contribution to the sum rules is 
expressed via the perturbative spectral density for which 
we use the corrected expression published as Eq. (1) in 
the Erratum to jtf], viz., 


Ppcr L t 0) (^) = ^r 


3xx 


1 + 


a s {p 2 

47T 


-C F 5 — 


- In 


(17) 

where C F = (N 2 - l)/2 N c = 4/3 for SU( 3) c . 3 The 
nonperturbative content of the sum rule is contained in 
the expression 


$ p(7r) (x, M 2 ) = =F$ 4 q(a;, M 2 ) + <S>q Aq (x , M 2 ) 
+$ v (x,M 2 )+<S> G {x,M 2 ), (18) 

where M 2 denotes the Borel parameter and sq marks in 
each case the threshold value which separates the lowest 
resonance state from higher states. To saturate the sum 
rules for the first N = 10 moments of T )(x) we 
use so ~ 2.25 GeV 2 . The various contributions, forming 
^ > p( 7 t); pertain to the following terms: (i) 4>4 q (four-quark 
condensate); (ii) <&qAq (quark-gluon-antiquark conden¬ 
sate) (iii) 4>v (vector quark condensate); (iv) $ G (gluon 
condensate). Their explicit expressions can be found in 
Appendix A of Ref. [481 ] . The basic assumption here 
is that higher-order correlations are less important than 
two-particle correlations (vacuum-dominance hypothesis 
f59i]). This assumption is employed in order to reduce 
the four-quark condensate to the product of two-quark 
condensates ignoring corresponding uncertainties. One 
notes that the four-quark contribution enters the sum 
rule for p || in (1181) with the opposite sign with respect 
to 7r. As a result, it reduces the relative weight of this 
condensate in the sum rule entailing smaller values of the 
DA moments - - in contrast to the pion case. In fact, we 
found in Q (see Table 1 there) the following relation: 

(£ 27V )tt > (£ 2W )ph > (£ 2Ar )asy with IV = 1, 2, 3. 

Similarly to the pion DA, the DA of the longitudinally 
polarized p is defined by the matrix element 

(0|J(z)7 M u(0)|p(p,A)| z2=o =/Jpp f dxe lx{z ' p) 

J 0 

x^pO^M 2 ), (19) 

while the definition of the transversal (_L) p DA reads 

(0\d(z)a^u(0)\p(p, A)| 22=q = ifp{e { p ) p v - e^W) 

x [ dxe’ lx( ' z ' p) ipj(x,p 2 ) , (20) 

Jo 

where we have again employed the gauge A+ = 0. The 
p _DA will not be considered in this work. 


3 This expression coincides with the 0(a s ) radiative correction 

entering the perturbative contribution to the sum rule considered 


G3- 


To construct <p\{x, p 2 ), we compute the moments 
(C Ar )py up to N = 10 from the QCD sum rule (fTTH) and 
determine from them the corresponding conformal coef¬ 
ficients a n with N = 0,2,..., 10. Because the moments 
with N > 6 turn out to have values close to the asymp¬ 
totic ones, we can safely set the associated conformal co¬ 
efficients equal to zero and express ipp(x,p 2 ) in the form 
of Eq. ([9]). In addition, we use the Gegenbauer-a expan¬ 
sion, given by Eq. m, via the parameters (a_, a% ). The 
accessible regions of these parameters for the longitudinal 
p DA, determined within QCD sum rules with nonlocal 
condensates, are displayed as shaded blue areas closer to 
the y axis in Fig. [2] The left panel displays the results 
in the ( 02 , 04 ) plane, while the right panel contains the 
analogous results in the (a_, 02 ) plane. We have also 
depicted in both panels the parameter regions referring 
to the pion case for A 2 = 0.4 GeV 2 (larger green slanted 
rectangles further to the right) and for A 2 = 0.45 GeV 2 
(transparent rectangles within dashed boundaries). The 
graphics in Fig. [2] include the platykurtic regimes of both 
DAs, determined in our present NLC-based analysis to 
be outlined later. Those areas referring to p\\ have red 
color and are located closer to the ordinate (left panel) 
and further to the top (right panel). The analogous ar¬ 
eas for the pion appear in light green color and are ad¬ 
jacent to the previous ones. In addition, we incorporate 
several other 7 r and p\\ DAs with individual designations 
explained in the figure caption for the readers’s conve¬ 
nience. The uncertainties of the presented pion and rho- 
meson DAs, obtained with our NLC QCD SR approach, 
include only those stemming from the SRs themselves 
and are related to the variation of the Borel parameter 
within the stability window of the SRs [3] . Experimental 
uncertainties in the input physical parameters have little 
influence on the results and have been ignored. 

The values of all parameters of the displayed models 
are listed in Table U at the reference scale p 2 = 4 GeV 2 . 
This scale is employed in lattice calculations because it 
naturally arises by the matching of the bare (lattice) op¬ 
erators at /Zq = 1 /a 2 (a being the lattice spacing) to 
those in the MS scheme in continuum QCD. It is also 
used in various works based on the DSE approach. In 
this work, we obtained our own results at the initial scale 
/i 2 > 1 GeV 2 and evolved them to this higher scale us¬ 
ing ERBL evolution at the next-to-leading order level. 
The conformal coefficients 02 and 04 , and the moments 
(£ 2 ) and (£ 4 ) along with the inverse moment (a: -1 ) in 
Table U| have for each DA their own original values at 
p 2 = 4 GeV 2 . In cases where the original scale was lower, 
we have evolved these quantities to the scale p 2 = 4 GeV 2 
using two-loop ERBL evolution. 

We have included in the table the AdS/QCD model of 
the pion DA derived in |43| within holographic QCD. 
This model reads <p T (x,p 2 ) = ( 4 /\/ 37 r)y/xE and ap¬ 
proaches at Q 2 —> 00 the asymptotic DA ip^ y (x) = 6xx, 
while it has a very different x behavior at finite Q 2 
EH 1>] • The conformal coefficients of this DA have 
been computed in the arXiv version of [20[ at the initial 
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TABLE I: Various parameters entering the DA Gegenbauer representations in m and G0 in the MS scheme. The kurtosis 
/? 2 , the second moment (£ 2 ), the fourth moment (£ 4 ), and the inverse moment (a; -1 ) of the n and the py DAs are also shown. 
Always the original functional forms of the DAs have been used with the same designations as in Fig. [2] The reference scale 
for all entries is p 2 = 4 GeV 2 either by construction or after NLO evolution. Only the central values of the normalization 
coefficients N a are displayed. Note that we ignore here and in the figures the numerically negligible effects on the asymptotic 
DA induced by NLO evolution [3ijl . 


Model DA 

ai 

a4 

N a 

a_ 

aj 

P 2 

<£ 2 > 


(z- 1 ) 

asy ♦ 

0 

0 

6 

1 

0 

2.14 

0.2 

0.086 

3 

P|| [18] ♦ 

0.056 

-0.055 

16.83 

1.63 

0.137 

1.93 

0.22 

0.093 

3.0 

pii [ 2 ] ★ 

0.032(46) 

-0.038(81) 

13.6 

1 50 +0 - 71 

± - ou — 1.50 

0 1 1 +0.08 
u - a± -1.14 

2.0(3) 

0.211(16) 

0.088(7) 

3.0(1) 

pjj 111 (here) © 

0.017(24) 

- 0 . 021 ( 11 ) 

10.0 

1 312 +0 ' 19 

0 071+ 0 - 016 
U.U/ l_o.oi5 

2.06(4) 

0.206(8) 

0.087(6) 

3.0(8) 

py [40] • 

-0.009 

-0.023 

11.81 

1.41 

0.063 

2.09 

0.197 

0.081 

2.92 

py [39] ▼ 

0.012 

-0.007 

7.18 

1.11 

0.034 

2.11 

0.204 

0.088 

2.98 

pii m ■ 

0.10 

0 

6 

1 

0.10 

1.98 

0.234 

0.109 

3.30 

pii m 0 

0.092 

0.031 

3.37 

0.66 

0 

2.05 

0.232 

0.110 

3.5 

PII [42] □ 

0.057 

0.002 

5.76 

0.975 

0.052 

2.05 

0.220 

0.099 

3.2 

7TBMS [ 6 ] X 

0 14Q + 0 - 052 
lJ.14ty_o.o43 

—0.096±g;gH 

20.49 

L - 1 °-0.45 

PI 91 7+0.048 
U.^i l —0.066 

1 74+0- 16 
*^-0.14 

0.248/^015 

0.108/°;°3 

Q i«+0.°9 
°- lu -0.09 

TTpk [5] 'fr 

0 n^ 7+ 0 - 024 
v.yjoi _o.oi9 

-0.013™ 

7.78 

1 1 a+0.24 
1-10-0.26 

0.086i°;°26 

9 no+0-02 
z.u^_o 03 

0 . 220 /“;°°® 

0.098/°;°o5 

0 1 q+0.14 
o.iO-o.io 

7TDSE-DB [9] A 

0.149 

0.076 

1.81 

0.31 

- 0.12 

2.0 

0.251 

0.128 

4.6 

ttdse-rl [9] V 

0.233 

0.112 

1.74 

0.29 

0.0029 

1.9 

0.280 

0.151 

5.5 

TAdS/QCD [43] A 0.107 

0.038 

2.55 

0.50 

0 

2.03 

0.237 

0.114 

4.0 a 

7TCZ [ 1 ] 

0.42 

0 

6.0 

1.0 

0.42 

1.54 

0.343 

0.181 

4.25 

TTlat [49] 

0.1364(154)(145) 

- 

- 

- 

- 

- 

0.2361(41)(39) 

- 

- 


“This value was obtained using only the first three terms of the 
conformal expansion in Eq. and is therefore not a precise esti¬ 
mate. 


scale p 2 » 1 GeV 2 to obtain (/^AdS/QCD = 0 250 and 
(£ 4 }AdS/QCD _ q 225 from which the first two conformal 
coefficients a ^ dS /QCD = 7/48 and a ^ dS /Q CD = 11/192 
were determined. Using NLO scaling relations, these co¬ 
efficients have been evolved to the reference scale p , 2 = 
4 GeV 2 and are given in Table Q] The value of the sec¬ 
ond moment was later computed in [50] for po = 1 GeV 
and Aqcd = 0.225 GeV using LO evolution and found 
to be almost the same, notably, = 0.24. 

Note that a broad pion DA ~ \fxi, was considered be¬ 
fore in ii using QCD sum rules with NLCs. An 
even broader DA was discussed earlier in !51], giving 
tp w (x) ~ (a;x)[° 1 _ 0 - 2 l. To complete this discussion, we 
mention that the coefficient a 2 (and other parameters) 
for the pion DA proposed by Chernyak and Zhitnitsky 
fl] has been included in Table Q] using NLO evolution to 
the scale p 2 = 4 GeV 2 , but we have not displayed it in 
Fig. [2] because its value is outside the range of the graphs. 

The most important observations from Fig. [2] are the 
following: (i) There is no overlap of the ( 02 , 04 ), or, 
equivalently, (a-, a?!), regions for the pion and the py 
meson computed with NLC sum rules in [@| and @]. (ii) 
While the platykurtic regime for the py DA is entirely en¬ 
closed by the region determined with A 2 = 0.4 GeV 2 , the 


platykurtic regime for the pion DA appears as an exclave. 
This is, because in order to obtain a platykurtic pion DA 
one has to use the slightly larger value A 2 = 0.45 GeV 2 . 
(iii) The DSE-based DAs for the pion and the py meson 
are in both panels far away from our NLC estimates. But 
keep in mind that the locations of the DSE DAs in the 
( 02 , 04 ) plane are only indicative because, as we have al¬ 
ready mentioned, their parametrization by means of Eq. 
© is a very crude approximation, (iv) We observe in 
the ( 02 , 04 ) plane an intriguing alignment of the unimodal 
DSE-based DAs ( 7 r and py) along an upward pointing “di¬ 
agonal” and another branch of bimodal NLC-based DAs 
steered downwards along 02 + 04 « const. This second 
“orbit” of DAs roughly follows the values of (a: _1 ) w /3 — 1 
evolved to the scale 4 GeV 2 @. Crucially, there is a small 
region where both branches overlap allowing the combi¬ 
nation of endpoint suppression with unimodality — the 
chimera regime of our interest, (v) On the other hand, 
all DAs, blended together according to Eq. (E6l) . lie on 
the straight line with 02 = 0.15 and stray away from the 
NLC ( 02 , 04 ) region with decreasing values of the mixing 
parameter a. Ultimately, i.e., for Q 2 —i 00 , all DAs with 
the correct x asymptotics will evolve either along the up¬ 
per “diagonal” (if they are unimodal) or along the lower 
one (if they are bimodal) toward the asymptotic DA (♦) 
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with 70 = 0, as predicted by perturbative QCD. Issues 
(iv) and (v) will be addressed later in full detail. 


V. MESON DAs AS PATTERNS OF 
SYNCHRONIZATION 

The meson DA at a fixed scale p 2 is a distribution 
of x values in the interval [0,1]. While a unimodal DA 
profile may seem more “natural” in appearance for the 
ground state of the pion, the interpretation of a bimodal 
structure causes discomfort and gives rise to debates. So 
there is a desire for an explanation and rationalization of 
this issue. 

Recently, it was argued by one of us 0] that in order 
to better understand the patterns of the pion and other 
meson DAs it is useful to develop some ideas which are 
drawn from the subject of synchronization of nonlinear 
oscillators in the theory of complex systems — natural 
and engineered (see {HU, 53] for reviews). The nub of 
the idea, as Stefanis put it in Ref. 3, is to represent 
the x values, accessible to the meson DAs in the interval 
[ 0 , 1 ], in terms of the phases of a large number (N —> oo) 
of interacting oscillators. The dynamics of this kind of 
systems is describable in terms of the Kuramoto model 
[54| and its descendants, albeit its specifics is not rele¬ 
vant for the present analysis. What is more important is 
that the synchronization of the oscillator phases, alias the 
longitudinal momentum fractions carried by the valence 
quark vs. that of the antiquark, gives rise to the forma¬ 
tion of particular patterns of the x distribution. These 
patterns emerge from the “organization” of the phase 
spectrum (i.e, the x distribution) and reflect the specific 
approach used to describe the partonic interactions in the 
pion bound state described by the DA. In other words, 
each particular DA profile is latent in the underlying the¬ 
oretical method and pertains to a patterned arrangement 
of synchronized coupled oscillators in the Kuramoto con¬ 
text. 

These methods can be QCD sum rules with nonlo¬ 
cal condensates, as employed in this work and in [®, 1 } 
(see also [ 22 ], [23) . local condensates [U■ lightcone sum 
rules (LCSR)s [23, [55], instanton models 4Cj, approaches 
based on Dyson-Schwinger-equations [l(| (reviewed in 
56| ), light-front models, e.g., [H, [57f, holographic QCD 
— see [HI for a recent review, etc. Examples of 7 r 
and p|| DA profiles are depicted in Fig.Q] Thus, the Sync 
concept provides a universal canvas to study the charac¬ 
teristics of very different meson DAs without taking re¬ 
course to a specific calculational scheme. In particular, 
it puts a theoretical basis beneath the interpretation of 
the bimodality of meson DAs, as we will show next. 

Indeed, it was pointed out in 3 that at scales of the 
order of /i ~ 1 — 2 GeV, nonlocal condensates, which are 
used to parameterize the vacuum nonlocality in terms of 
a nonvanishing quark virtuality A 2 cause the distribution 
over x to flock into two clusters, giving rise — within a 
broad range of uncertainties in the midregion of x — to 


bimodal DAs for the pion [b] and the p\\ meson Q. The 
corresponding families of DAs are shown in the form of 
the larger shaded bands in Fig. [T] The upper panel refers 
to the 7r and the lower one to the p\\ meson, both at the 
scale p 2 = 4 GeV 2 after NLO evolution from the initial 
scale Pq ~ 1 GeV 2 . 

The bimodality strength of the BMS-type of DAs is 
controlled by the nonlocality parameter A 2 , which en¬ 
dows vacuum fluctuations with a characteristic corre¬ 
lation length ~ 1/A q . Lower values of A 2 tend to in¬ 
crease the bimodality character of the DA and reduce 
the value of '■p 7l {x = 1/2), while larger values enhance the 
midregion of x driving ip n (x) closer to a unimodal distri¬ 
butional shape. This behavior is deeply rooted in the 
combined effect of the perturbative part and the power- 
behaved terms in the QCD sum rule for the moments 
(£ N ) considered in Q and in [§]. For A 2 = 0, one recov¬ 
ers the QCD sum rules of Chernyak-Zhitnitsky in Ref. 
[I] with an infinite correlation length of the vacuum fluc¬ 
tuations. The numerically most important term is the 
scalar-condensate contribution [l3j, encountered in (fl8l). 4 * 

(x;M 2 ;A) J^{6»(5 > A > x)x [x 

+ (A — x) In (x)] + (x —> x) +0(1 > A) 
x 9 (A > x > A) [A + (A — 2xx) 
x ln(A)]} , (21) 

where Ag = (8na s /81 )(qq) 2 , with the four-quark contri¬ 
bution being given by a s (qq) 2 = 1.83 x 10 -4 GeV 6 and 
(a s GG)/127r = 0.0012 GeV 4 Here, A = A 2 /2M 2 , 
A = 1 — A, and M 2 ss 1 GeV 2 is the Borel parameter 
M 2 6 [M r 2 in , M/ lax ]. The sum rule should not be sensi¬ 
tive to the choice of this parameter. The procedure for 
minimizing the dependence on M 2 has been described in 
[H| and is applied here. 

Larger values of A shift the balance in the sum rule in 
favor of the perturbative contribution which has a sin¬ 
gle mode at x = 1/2, thus entailing a reduction of the 
two peaks at xo = A and Xo = A until they ultimately 
collapse into a single more rounded peak at the center. 
But despite this shift towards the central region of x, the 
tails of the BMS DAs remain suppressed within only a 
small range of theoretical uncertainties, as one also ob¬ 
serves from Fig. [T] Moreover, it was shown in J48| that 
the endpoint behavior of the pion DA can be related to 
the “decay rate” of the correlation length of the scalar 
condensate. It is worth noting in this context that in the 
local version of the condensate model, i.e., for A 2 = 0, all 
nonperturbative contributions are concentrated just at 
the endpoints because 4>] l < Q al (x) = 9 [<5(x) + <5(x)] /(M 2 ) 2 . 

Viewed from the Kuramoto prism, the two peaks of 
the BMS DAs correspond to two distinct groups of rather 


4 The expression for the scalar quark-condensate contribution in 

Eq. (1211) pertains to a Gaussian model for the quark condensate. 






strongly synchronized oscillators with characteristic “fre¬ 
quencies” located in the lower and upper quartiles of the 
x distribution, respectively. On the other hand, the tails 
correspond to tiny cohorts of oscillators with natural “fre¬ 
quencies” close to the rare values x = 0 and x = 1 with 
almost nil phase-locking, while a partly synchronized ar¬ 
rangement of oscillators with values around x = 1/2 con¬ 
nects the two clusters across the dip in the central region. 

Consider now the implementation of DCSB to meson 
DAs. It was stated in |4| that the DCSB and the con¬ 
comitant mass generation of quarks and gluons within 
the DSE-based framework [§, |To| tend to enhance, both 
the central region of x values but also the tails of the me¬ 
son DAs down to the kinematic endpoints x = 0, 1, lead¬ 
ing to a homogenization of the x values of the valence 
qq pair and to broad unimodal DAs for all considered 
mesons GS H3, H2 • These DAs have downward con¬ 
cave profiles in the whole interval x £ [0,1]. There are 
two variants of pion DAs derived from the DSE-based 
approach [9] . They were computed via a large number of 
moments ((x — x) N ) ( N = 50) and were then expressed 
by means of Eq. (fT2l) using two different DSE truncations 
at the scale /r 2 = 4 GeV 2 . The associated values of the 
parameters (N a , a_, ) and other relevant metrics are 

provided in Table [U 

One 7r DA, dubbed DSE-RL, was obtained using 
the rainbow-ladder (RL) approximation of the Bethe- 
Salpeter kernel in the DSEs, while the other, termed 
DSE-DB, was derived with a DCSB-improved kernel 
(abbreviated by DB), which includes nonperturbative 
DCSB-generated effects that were not taken into account 
in the RL truncated version. Both DAs are much broader 
relative to (f^ sy , with the DSE-RL DA being flatter and 
broader than the DSE-DB DA. The profiles of these 
DAs are given in the upper panel of Fig. [T] DSE-RL 
dashed-dotted line; DSE-DB — dashed line. Also the pn 
DA obtained with the DSE computational scheme [lOj] , 
is a relatively broad everywhere downward concave curve 
(though less pronounced than both pion DSE DAs), as 
one can see from the lower panel of Fig. [T| (lower dashed 
line). The broadening of the DSE DAs is a direct conse¬ 
quence of the nonperturbative DCSB interactions which 
give rise to the dressed quark’s selfenergy — see [56| for a 
detailed review of the method and further explanations. 
Also note that the DSE DAs cross the NLC-based ones 
twice on each side of the mean (Fig. [TJ. Thus, these 
DAs show a pattern of higher-lower-higher on each side, 
related to their heavier tails. This behavior can be quan¬ 
tified by employing the kurtosis statistic, defined by 

E ( x - /x[y>]) 4 = /mM = (£ 4 )tt , 22 x 

{E{x-„[v]y? **&] ((Ox) 2 ’ 



mesons composed of light and heavy quarks like the kaon) 


7i M 


E(x - n[<p}) 3 

a 3 [tp\ 


(24) 


it describes the central tendency, variability, and shape of 
a distribution. In particular the kurtosis serves to mea¬ 
sure the peakedness in the central region of a distribution 
against the flatness of its tails. As one sees from Table 
[H the unimodal, downward concave DSE-based DAs can 
be ordered as follows: > @ 2 DB > P 2 RL ■ This result 

confirms a similar qualitative statement in M- 

Such broad DA morphologies describe a loosely syn¬ 
chronized assortment of oscillators spread over the entire 
range of their native “frequencies” in £ £ [0,1]. The en¬ 
hancement of the generic midregion, which corresponds 
— from a physical perspective — to “egalitarian” par- 
tonic configurations in the pion in which the valence 
quark and the valence antiquark carry comparable frac¬ 
tions of longitudinal momentum, may be welcome. But 
at the same time DSE DAs also overestimate the weight 
of “aristocratic” configurations with low probability oc¬ 
currence in which a single valence parton takes the lion’s 
share of the momentum with x —> 1 or x —> 1 to go 
far-off shell. Because these are more specific and rare 
configurations of the dispersion of the valence parton’s 
longitudinal momentum far away from the typical val¬ 
ues around the mean /r = 1/2, tail enhancement is at 
odds with our understanding of the QCD description of 
the pion bound state based on off-shell gluon exchanges 
(see [62], H and [50] for explanations) and, as we have 
seen, leads to more variation in the oscillator “frequen¬ 
cies” , thus entailing less synchronization. Note that the 
broader and flatter a unimodal downward concave x dis¬ 
tribution is, the closer to random the oscillator phases 
will be. 

Comparison in earlier works, e.g., sain, of predic¬ 
tions for the pion-photon transition form factor obtained 
within the LCSR framework with all existing data, in¬ 
dicates that strict QCD scaling behavior at high Q 2 is 
very sensitive to the endpoint-behavior of the pion DA. 
This behavior is intimately related to the inverse moment 
(x _1 ) w . Recalling Eq. m, this implies that the confor¬ 
mal coefficients have to balance each other in such a way 
so that the excess over (a: _1 )^ sy = 3 is not too large. Oth¬ 
erwise, the form factor would overshoot the data. And 
in fact, using the DSE DAs as nonperturbative input in 
a LCSR-framework, we have shown j|| that the predic¬ 
tions obtained herewith exceed the asymptotic limit even 
at the highest momenta around 40 GeV 2 probed in cur¬ 
rent experiments [i(| |4^ . This is also obvious from Table 
IU The inverse moment 


where 

r 1 i 

° 2 W\ = J o dx{x - fj,[ip]) 2 <p(x) = -(£ 2 ),r (23) 

is the variance of the distribution ip(x). Together with 
the skewness (vanishing here but being relevant for 


{x -% = i±^(l + a « + ...) , (25) 

a- 

obtained with both pion DSE DAs, has very large values 
to be compatible with the data. We mention in this con¬ 
text that a recent analysis jU of the meson structure in 







9 


lightfront holographic QCD, which employs a broad uni- 
modal 7r DA, finds predictions for the pion-photon transi¬ 
tion form factor that are in quite poor agreement with all 
experimental data (see Fig. 18 there). But notice that 
a subsequent DSE-based computation :66j of the TFF, 
in which the role of the inverse moment is not so crucial, 
finds good agreement with the CELLO, CLEO, and Belle 
data in the entire domain of spacelike momenta. 

In the Sync analogy, the NLC within our approach 
causes a generic clustering of the DA into two clusters 
liaised with massive endpoint suppression - within the 
range of intrinsic theoretical uncertainties — and only a 
moderate reduction of the DA in the central region, which 
is controlled by the strength of the nonlocality parame¬ 
ter as discussed above in connection with Eq. ET1) . Thus, 
configurations with a highly asymmetric dispersion of the 
valence parton’s longitudinal momentum are suppressed. 
It is highly unlikely that the pion can remain intact and 
rebound as a whole for such partonic configurations. In 
the Sync analogy they would play the role of very idiosyn¬ 
cratic oscillators incapable to synchronize, being either 
too slow (with “frequency” values close to 0) or too fast 
with “frequencies” tending to unity. In this sense, the 
NLC acts like a negative feedback opposing the excessive 
DCSB-induced enhancement in the endpoint regions by 
turning off the corresponding oscillators. A figurative ex¬ 
planation of these issues is provided in Fig. [3] in which 
the two main antithetic effects in the x behavior of the 
valence quark in the pion DA are illustrated. This is how 
Stefanis envisioned and mapped out the Sync properties 
of the BMS-like and DSE-like pion DAs in Q. So there 
are, it seems, two very distinct DA patterns with tell¬ 
tale signatures that include, but are not limited to, the 
pion case: One is related to NLCs, which encode par¬ 
ticle correlations in the range 1/A g ~ 0.3 fermi, while 
the other implements DCSB which causes the dynamical 
generation of quark masses and entails dressing of the 
quark propagator describing the confined quark in the 
pion. Both effects are manifestations of confinement and 
neither exists in isolation — see [67] for a recent review 
of strong-interaction dynamics. Unfortunately, they can¬ 
not be described simultaneously within a single analytic 
approach at present. Thus, for the time being, there are 
two physical paradigms with their own computational 
techniques, each applying only within its own sphere of 
acceptance and validity. 


VI. SYNTHETIC MESON DAs 

There are basically two options: (i) Either the x dis¬ 
tribution of the pion DA is described by a single DA over 
the whole range of values in the interval x € [0,1], or (ii) 
the “true” pion DA is rather a mixture of two different 
DAs, one better applicable to the central region and the 
other controlling the tails. The first option was discussed 
above and at length in the literature. Following the sec¬ 
ond scenario, Stefanis proposed in [3| a synthetic DA of 



FIG. 3: (color online). Illustrating the key antithetic effects in 
the pion distribution amplitude in the interval [0.5, 1] (mirror 
graph in the interval [0,0.5] not shown). The lower (green) 
line denotes the BMS DA, with the range of such DAs being 
indicated by the green vertical strip (cf. green shaded band 
in the upper panel of Fig. [I]). The upper (blue) line marks 
the DSE-DB DA in the central x region. The unimodal DSE- 
DB DA shows enhancement of the tails and the midregion 
around x = 0.5, while the BMS 7r DAs are characterized by 
suppression of the endpoint region x = 1. The correspond¬ 
ing main trends are indicated by vertical arrows: DCSB 
enhancement; NLC — suppression. 


the form 

rf ue (x : a) = (up® MS (x) + (1 - «)^ SE (:r) , (26) 

where a is a mixing parameter with values within the in¬ 
terval [0,1]. Mixtures of the form of Eq. (l26l) are quite 
common in statistics when a single distribution, like the 
Gaussian distribution function, the Poisson distribution, 
etc., has to be combined with a distribution with a dif¬ 
ferent type of mathematical behavior in the tails, e.g., 
the generalized Pareto distribution. In fact, hybrid-like 
DAs have been constructed and profoundly studied by 
Bergmann and Stefanis long ago for the nucleon [68l — 
Iz3 and also for the A + (1232) resonance (termed 
“heterotic” DAs), although they were motivated by other 
concerns. A comprehensive and detailed review of such 
baryon DAs has been given in Q. In the present case, the 
coexistence of distinct domains of oscillators, some coher¬ 
ent and phase-locked (BMS peaks), and others which de¬ 
scribe unsynchronized oscillators (heavy tails of the DSE 
DA), would give rise to a so-called chimera state jut [73} ■ 
It was argued in Q that for values of the mixing coeffi¬ 
cient a close to 1, the synthetic DA would still belong to 
the family of pion BMS DAs shown in terms of the wide 
shaded band in Fig. [T] More generally, the synthesized 
DA is supposed to encapsulate both manifestations of 
QCD confinement, reflecting the perpetual balance which 
arises from the appropriate combination of the two ba¬ 
sic effects, one associated with NLC formation, via QCD 
sum rules (BMS-DA @), the other being the result of 
DCSB, expressed in terms of a DA computed within the 
DSE-based approach Q. A proper combination has to 
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FIG. 4: (color online). Synthetic pion DAs obtained with Eq. 
(1261) at the scale p 2 = 4 GeV 2 for various values of the mixing 
parameter a € [0,1], The shorttailed platykurtic (pk) DA is 
shown by the upper red solid line for comparison; it does not 
belong to this class of DAs (see text). 


balance the enhancement impact of DCSB against the 
suppression due to NLC, as exposed in Fig. [3] 

Under ERBL evolution the synthetic DA would de¬ 
velop at Q 2 —> oo to the asymptotic DA which repre¬ 
sents the most synchronized configuration of the valence 
qq pair within the pion being still a bound state [4| after 
all quark-gluon interactions have died out. This is also 
evident from Table [T] from which we see that the asymp¬ 
totic DA has the largest kurtosis, i.e., </? asy is the most 
leptokurtic meson DA. In this paper, we work out Eq. 
m in more certain terms and exploit the whole range 
of possible values of a. From the synchronization point 
of view, a synthetic DA represents an attempt of combin¬ 
ing ensembles of synchronized and unsynchronized (but 
otherwise identical) oscillators in order to enhance or sup¬ 
press particular frequency values amounting to a chimera 
state. The question is whether the simple one-parametric 
design of Eq. (E6l) is indeed capable of providing the de¬ 
sired properties addressed above in the adjunct discus¬ 
sion of Fig. O 

FigureS] shows various samples of synthesized DAs ob¬ 
tained with different values of the mixing parameter be¬ 
tween 0 and 1 . As one clearly sees from these plots, when 
a is larger than, say, 0.75, the main characteristics of the 
BMS DAs persist. This is also obvious from the loca¬ 
tion of this DA, denoted by the symbol (©), in the plane 
( 02 , 04 ) in the left panel of Fig. [2] However, for small a 
values close to 0.25 and below, the obtained DA profiles 
show no endpoint suppression, as desired, but exhibit in¬ 
stead tail enhancement, inherited to them by the DSE 
DA, while the profile is still bimodal, see ® (o = 0.25) 
and ® (a = 0.5) in the same figure. This means that 
the bimodality of the DA prevails from large down to 
quite small values of the mixing parameter while at the 
same time for these small values of a the tails of the DA 
get strongly enhanced. Hence, Eq. (Edrl) cannot supply a 
pion DA which combines unimodality in the central re¬ 
gion with suppression of the tails, despite the fact that it 


generates DAs with the same value of 02 , see, Fig. [2J left 
panel. Nor can this be realized via ERBL evolution . 5 In 
order to embody endpoint suppression of the pion DA, 
one has to build it in right from the start, resorting to 
QCD sum rules with NLCs and looking for DAs which 
would mimic the characteristics of the DSE DAs in the 
central region, while preserving suppression of the tails in 
compliance with Fig. [3] To achieve this goal we have to 
keep in mind that if mass is moved from the shoulders to 
the center of a distribution, then one has to compensate 
the accompanying movement of mass to the tails, leaving 
the variance almost unchanged but increasing the kurto¬ 
sis. 

The existence of such a chimera DA for the pion, which 
binds these diverse aspects of coherence and incoherence 
into a single DA, was first discussed in Q and predictions 
for the pion-photon transition form factor Q 2 F J 77r ( Q 2 ) 
were presented which are fully compliant with all exper¬ 
imental data compatible with QCD scaling. The overall 
quality of these predictions resembles that of the BMS- 
type DAs 0, @]. This can be traced back to the fact that 
they both lead to an inverse moment with just the appro¬ 
priate size in order to agree with the Q 2 F 1 77r ( Q 2 ) data. 
In fact, a brand-new simultaneous fit to the CLEO (45| 
and Belle [13 data in [65] favors a profile of the pion DA 
which is very close to the platykurtic one (see Fig. Q] and 
TablelU). On the other hand, DAs with downward concave 
shapes in the entire interval x £ [ 0 , 1 ], will tend to over¬ 
estimate most data of the Belle Collaboration [13], the 
reason being that they giv e rise to large inverse-moment 
values (Table QJ and Hj5] (see also [ 6 J] ). [The opposite 
behavior was found in | 66 [, as already mentioned.] 

In this work, we determine a whole domain of such 
chimera 7 r DAs employing our NLC technology and al¬ 
lowing for a slightly larger value of the quark virtuality, 
viz., A 2 = 0.45 GeV 2 . The core area of these DAs is il¬ 
lustrated in the upper panel of Fig. [I] in the form of a 
narrow strip in red color. The corresponding parameters 
and the range of theoretical uncertainties for both used 
Gegenbauer parametrizations are given at the reference 
scale n 2 = 4 GeV 2 in Table Q] This table also includes 
the brand-new lattice estimates for the second moment 
(£ 2 ) from [HI at the same scale using the MS scheme. 

Note that this 02 value was not calculated from the 
second moment (£ 2 ) via Eq. (0 at finite lattice spac¬ 
ing, taking subsequently the continuum limit. Instead, 
the value of 02 was calculated directly on the lattice 
and was then extrapolated to the continuum limit via 
(£ 2 )!^o =*■ afWo => a“ S |a=o- This implies that Eq. © 


5 Strictly speaking, evolution in NLO gives a logarithmic modifi¬ 
cation ~ (a s /(47r))CF ln 2 [(l —x)/x\ 1331 which obviously affects 
the endpoint behavior of the meson DA — independently of its 
shape — albeit endpoint-enhanced DAs receive larger NLO cor¬ 
rections [53 . However, given that at scales /z > 1 GeV 2 the 
running coupling is already sufficiently small, this effect can be 
safely neglected, see 0. 













11 


is broken by lattice artifacts. A small variation in the lat¬ 
tice spacing around 6 % may result in an increase of a 2 of 
the order of 25—30% [49j. The final result at p 2 = 4 GeV 2 
reads a^ s | a= o = 0.1364(154)(145), while the reported 
value of the second moment is (C 2 )^* = 0.2361(41)(39). 
The first error is statistical and originates from the chi¬ 
ral expansion, whereas the second one pertains to the 
uncertainties of the renormalization factors. It agrees 
within errors with (£ 2 ) BMS = 0.248{g;|g| [H] and also 
with (£ 2 )£ k = 0 . 220 {q ;| 2 3 (see TableO}. In contrast, while 
the lattice estimate a 2 at agrees with the BMS coefficient 
af MS , determined in the year 2001 0 (see Tabled]), it 
turns out to be larger than a 2 k . But one should be cau¬ 
tious. Extracting a 2 via (C 2 )^f 0 (£ 2 )^o => a^ s | a=0 , 

one would obtain a 2 at = 0.105, which is indeed compati¬ 
ble with the range of the platykurtic a 2 values. Thus, one 
cannot exclude the influence of significant discretization 
effects that would require simulations at smaller lattice 
spacings of the order of a ~ 0.04 fm El] , 

Be that as it may, one should recall that the second 
moment (t; 2 )^ is related to the variance of the DA given 
by Eq. (l23l) . This statistic is not sufficient to draw any 
conclusions about the shape of the distribution in the 
central region. Indeed, as one observes from Table dl 
the unimodal DSE-DB pion DA yields a conformal co¬ 
efficient a^ SE-DB = 0.149 which fully agrees with the 
new lattice result but also with af MS . We note that this 
is valid for the second moment as well, which has the 
value (£ 2 ) ESE-DB = 0.250, (cf. (0) and thus almost co¬ 
incides with the second moment of the pion BMS DA 
given above, being also close to (£ 2 )£ k . On the other 
hand, the fourth moment (£ 4 ) and the conformal coeffi¬ 
cient 04 of the DSE and the BMS (pk) DAs are different 
in value and sign, respectively — see Table H] What is far 
more significant is the fact that, as it is evident from the 
left panel of Fig. [2] all DAs lying on the straight vertical 
line at o 2 =0.15 agree equally well with the new lattice 
estimate for o 2 . This makes it apparent that a single 
lattice constraint cannot fix the profile of the pion DA 
uniquely, however precise it may be. 


The chimera DAs have shorttailed platykurtic profiles 
and overlap with the DSE-DB DA in the midregion of x 
but descend at the endpoints at low angle to zero, similar 
to a typical BMS DA. As one observes from Fig. [2] (both 
panels), there is an imbrication of the platykurtic regimes 
(small rectangles in light-green color) with the domains 
of the bimodal pion DAs obtained with NLC sum rules 
for the quark virtuality A 2 = 0.45 GeV 2 (transparent 
rectangles bounded by a dashed line). For this value the 
conformal coefficients for <^ BMS (;r) at p 2 = 4 GeV 2 read 
a 2 = 0.12 and 04 = —0.06 while the inverse moment 

is (ar 1 )™ 3 ^*” 0 ' 45 GeV2) = 3.18, a value which agrees 
well with (x _ 1 )P k = 3.13 in Table U] The prediction for 


Q 1 F 1 77 r °(Q 2 ) obtained with the shorttailed platykur¬ 
tic 7 r DA appears in line with all data of the Belle ^7| 
and the BaBar Collaboration [46j compatible with QCD 
scaling 30- These unique features of the pk pion DA 


look indeed very attractive. But is it more than mere 
coincidence or can it provide a general mode of accessing 
meson DAs and offer a deeper perspective on meson DAs 
in general? 

To this end, we turn our attention to the p meson case 
and attempt to determine a platykurtic regime for the 

a I DA using as a selector the behavior illustrated in Fig. 

Evaluating the sum rule in Eq. (USD, we compute the 
reliability range of the conformal coefficients up to the 
order N = 10 by first determining the central moments 
(£ N )p\\ °f the same order. Their values at the initial scale 
p 2 > 1 GeV 2 can be found in 0. Also the correspond¬ 
ing values of the p|| meson are given there together with 
the conformal coefficients. We will not repeat these de¬ 
tails here. We concentrate instead on our primary goal 
to extract a platykurtic domain of these parameters. It 
turns out that this is possible even for the somewhat 
smaller value of the quark virtuality A 2 = 0.4 GeV 2 used 
originally for the extraction of the pion DA in 0. The 
extracted domains are shown in Fig. [2] in the form of 
the red slanted rectangles surrounded by the larger blue 
bands of coefficient values computed with the NLC sum 
rules. The left panel displays the results for the first two 
conformal coefficients a 2 and < 24 , whereas the right panel 
provides the areas of the coefficients a_ and a% ■ In both 
graphics the platykurtic model for the p\\ DA is denoted 
by the symbol ©. The values of all these parameters, ac¬ 
companied by their intrinsic errors, are compiled in Table 
HI while the platykurtic p\\ DA profiles are displayed in 
the form of a narrow red strip in the lower panel of Fig. 
[U For the sake of direct comparison with the DSE re¬ 
sults, all graphics and the values in the table are given 
at the reference scale p 2 = 4 GeV 2 after two-loop evo¬ 
lution. One immediately observes from this figure that, 
similar to the pion case, the platykurtic p\\ DA has a sin¬ 
gle rounded central peak bearing endpoint suppression 
relative to the p\\ DA obtained within the DSE-based ap¬ 
proach 0 - In comparison to the platykurtic DA of the 
pion, it features a slightly narrower profile with the kur- 

^pk k 

tosis value /3 2 11 > /3J . Tangible consequences of the 

platykurtic p» DA will be studied elsewhere. 


VII. CONCLUSIONS 


We have performed an intensive study of the pion and 
the p\\ DAs within QCD, fortified with the knowledge 
of synchronization concepts used in the description of 
complex systems. These concepts provide a unifying ra¬ 
tionale of how the various DA profiles emerge instead 
of asking why they should have a particular shape, thus 
avoiding descriptive comparisons of DAs obtained with 
unrelated theoretical frameworks. Furthermore, guided 
by these concepts, we have used controlled theory tools 
to obtain a new kind of chimera DAs for the pion and the 
PH meson using QCD sum rules with NLCs. These DAs 
are capable of mingling in situ the best of both worlds - 
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endpoint suppression via NLC and unimodality due to 
DCSB, giving rise to shorttailed platykurtic profiles and 
realizing the scenario illustrated in Fig. O while preserv¬ 
ing the asymptotic x behavior predicted by perturbative 
QCD. In the Sync picture, they correspond to a vast num¬ 
ber of phase-locked oscillators between the lower and the 
upper quartile of the x distribution, whereas oscillators 
with extremely high or low “frequencies”, located close 
to the tails x = 0,1, are in limbo. While the character¬ 
istics of these new DAs in the central x region resemble 
the gross behavior of DSE-based DAs, their suppressed 
tails are following the same trend as the BMS DAs. In 
mathematical terms, the BMS-like DAs and the platykur¬ 
tic ones are very different as regards their profiles (Fig. 
[T]) and Gegenbauer coefficients (Fig. [2] and Table Q}. But 
from the NLC point of view, the bimodal BMS n DAs and 
the bimodal p|| DA of which in the Sync analogy cor¬ 
respond to two clusters of synchronized oscillators, are on 
the same theoretical footing as the unimodal shorttailed 
platykurtic DAs for these mesons, which unite the phase- 
locked oscillators in a single group. Moreover, in the pion 
case they yield coinciding predictions for the pion-photon 
transition form factor which agree well with all available 
experimental data compatible with QCD scaling above 
~ 9 GeV 1 2 3 4 5 6 7 8 9 . Given all these results, we don’t want to 
stretch the importance of unimodality too far. 

Too broad DAs with downward concave profiles en¬ 
compassing the tails, as those derived for mesons with 
the aid of DSEs [§, [Tcj , imply that there is no particu¬ 
lar x value standing out because even the remote regions 
close to the endpoints x = 0,1 have a significant weight 
almost comparable to that of the central region — espe¬ 
cially the DSE-RL pion DA. The extreme case of a flat- 
top DA, like <z^~ top (x) = r(2(a + l))[r 2 (o!+l)] _1 (a;x)“ 
with a = 0.1 jig, translates into a population of oscilla¬ 
tors with a very strong variation of native “frequencies” 
so that these oscillators can hardly synchronize and as 


a result phase locking diminishes. Physically, this kind 
of x distribution comprises extremely asymmetric par- 
tonic configurations that can spoil scale locality and thus 
collinear factorization. On the experimental side, flat-top 
DAs yield predictions for the scaled pion-photon transi¬ 
tion form factor which have a tendency to increase with 
Q 2 — at least in the domain of currently accessible mo¬ 
mentum values in the range 10-40 GeV 2 where one would 
expect scaling to be visible j28(. The high-Q 2 data of 
the BaBar Collaboration [46[ indicate such a trend, but 
are not supported by the Belle data (T7j in the same re¬ 
gion. The next-generation experiments to measure the 
pion-photon transition form factor with the Belle II de¬ 
tector at the upgraded KEKB accelerator (SuperKEKB) 
in Japan and more precise data on the electromagnetic 
pion form factor expected at the Jefferson Laboratory 
(JLab) after its upgrade will provide extraordinary tools 
to test our predictions and assertions. 
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